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ABSTRACT 

cn 

^ ' Context. The possible existence of a continuum encompassing the diversity of explosive stellar deaths, ranging from ordinary super- 

' novae (SNe; lacking any sign of a relativistic outflow) to relativistic hypernovae associated with energetic long duration gamma-ray 

, bursts (GRBs), is under intense debate. In this context, the supernova SN 2008D associated with the X-ray transient (XT) 080109 could 

. represent a paradigmatic case, since it might exemplify a potential borderline transition event. Optical polarimetric studies could con- 

' tribute to shed light on the dilferent interpretations given in the literature for this supernova (hereafter, SN 2008D/XT 080109). 

Aims. The main aim is to infer geometric information of SN 2008D/XT 080109 through the study of the evolution of its linear optical 
polarization. We also report the polarization evolution of SN 2007uy, and discuss the properties of the host galaxy interstellar medium 
(ISM) towards the XT. The final goal is to compare the polarization properties, and therefore the geometries, of both SNe. 
Methods. We present a V-band linear polarization monitoring campaign carried out for SN 2008D/XT 080109 and SN 2007uy, which 
shone for weeks contemporaneously in NGC 2770. This fortunate coincidence brought us the opportunity to observe both objects 
simultaneously, and most importantly, with identical instrumental setups. The observations span 74.9 days, starting 3.6 days after 
the XT and are distributed in 1 1 visits. In addition we performed observations in the millimetre (mm) range in order to identify the 
■ dominant origin of the observed polarization. 

$H ' Results. We report positive linear polarization detections at several epochs for SN2008D/XT 080109 at a level of ~1%. For SN2007uy 

the measured polarization is around ~1.5%. In both cases the observed linear polarization seems dominated by the host galaxy in- 
terstellar polarization (HGIP), especially for the case of SN 2007uy. SN 2007uy shows Stokes parametres consistent with no time 
evolution, which could be described by the HGIP plus a constant eccentricity expansion on the sky plane. Over the course of our 
observations of SN2007uy, we find that its total polarization signal does not change by more than 0.29% with a 90% confidence 
interval. Despite the dominant HGIP, a statistical analysis of the distribution of the SN 2008D/XT 080109 Stokes parametres suggests 
that it could show a possible intrinsic variable polarization component. Moreover, assuming the polarization signal from SN 2007uy 
is constant, we find that the temporal evolution of the intrinsic SN 2008D/XT 080109 polarization could be explained by an aspherical 
axisymmetric expansion with variable eccentricity, although other more complex geometric scenarios are also compatible with the 
data. We come to the same result even if we make no assumption on the SN 2007uy Stokes parametres, although at a lower signifi- 
cance level. 

Conclusions. We conclude that the data seem to suggest a potential symmetry axis for SN2008D/XT 080109, which is reinforced 
when SN 2007uy is assumed to have constant Stokes parametres and used as reference star. We suggest that at least the projected, if 
not the intrinsic, geometry of SN 2008D/XT 080109 and SN 2007uy could be difi'erent. 

Key words. Stars: supernovae: general - Stars: supernovae: individual: SN2008D/XT 080109 ~ Stars: supernovae: individual: 
SN2007uy - Gamma rays: bursts ~ Techniques: polarimetry 



1. Introduction 

The X-ray transient (XT) 080109 was serendipitously discov- 
ered on 2008 January 9.56446 UT by Swift as a bright X- 



ray outburst during a scheduled observation of SN2007uy 
(Berger & Soderberg 2008; Kong & Maccarone 2008). Optical 
follow-up observations revealed the presence of a supernova 
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(SN) associated with the XT (Deng & Zhu 2008, hereafter 
named SN2008D/XT080109). SN2008D/XT080109 under- 
went a transition 1-2 weeks after the XT from Type Ic to 
lb (Malesani et al. 2009; Blondin et al. 2008; Valenti et al. 
2008; Modjaz et al. 2009), resembling the spectral evolu- 
tion seen in SN2005bf (Anupama et al. 2005; Modjaz et al. 
2005; Tominaga et al. 2005; Wang & Baade 2005). Eventually, 
SN2008D/XT080109 was spectroscopically classified as a He- 
rich lb supernova (Modjaz et al. 2009). 

The origin of this XT is still unclear. Two main scenarios 
have been proposed to explain its nature; the stellar shock break- 
out (Soderberg et al. 2008; Chevalier & Fransson 2008) and the 
X-ray flash framework (Li 2008; Xu et al. 2008; Mazzali et al. 

2008) . In the first case the early emission would be entirely ex- 
plained by the supernova, whereas in the second case a relativis- 
tic outflow from a central engine would be required. Radio data 
seem to indicate no superluminal motions, arguing against the 
presence of a long-lived relativistic outflow (Bietenhloz et al. 

2009) , therefore supporting the first framework. However, dif- 
ferent spectral analyses of the X-ray data did not yield conclu- 
sive results on the possible presence of a thermal component 
expected from the shock breakout wave (Chevalier & Fransson 
2008; Li 2008; Soderberg et al. 2008; Modjaz et al. 2009). 
Witiiin the second scenario SN2008D/XT 080109 was proposed 
as a transition between the most energetic hypernovae (linked to 
canonical long duration gamma-ray bursts) and standard core- 
collapse SNe (Mazzali et al. 2008). In fact it was suggested that 
SN2008D/XT080109 is a side-viewed, bipolar explosion with 
a viewing angle of > 50° from the polar direction (Tanaka et 
al. 2009a). It should be noted that spectroscopic observations of 
SN 2008D/XT 080109 suggested an aspherical expansion geom- 
etry (Modjaz et al. 2009; Tanaka et al. 2009a). 

SN2007uy was discovered on 2007 Dec. 31.669 UT 
(Nakano et al. 2008). It was classified as an ordinary lb SN 
(with no hint of any superluminal outflow) before the light curve 
maximum (Blondin & Calkins 2008), showing spectral similar- 
ities wifli SN2004gq (Pugh et al. 2004; Modjaz & Falco 2005). 
In the optical botii SN 2008D/XT080109 and SN 2007uy shone 
contemporaneously for weeks in the nearby (z = 0.007) spiral 
galaxy NGC 2770. This fortunate coincidence brought us the op- 
portunity to observe both objects simultaneously, and most im- 
portantly, with the same instrumental setup. 

Recently Maund et al. (2009) have reported a spectropolari- 
metric sfiidy of SN2008D/XT 080109 based on data acquired 
on 2008 January 31.22 and February 15.18 UT. One of the main 
advantages of spectropolarimetry of SNe with respect to broad- 
band polarimetry is its ability to infer geometric and dynamical 
information for the difl'erent chemical constituents of the explo- 
sion. Broad-band polarimetric observations construct a rougher 
picture of the stellar death, but require lower signal-to-noise ra- 
tios than spectropolarimetry. Hence broad-band polarimetric ob- 
servations can be extended to objects at higher redshifts or/and 
they allow to enhance the polarimetric coverage and sampling of 
the light curve, especially at epochs far from the maximum when 
the SN is dimmer. 

The dominant intrinsic polarization in SNe is thought to orig- 
inate by Thompson photon scattering through an aspherical pho- 
tospheric expansion (Hoflich 1991). If the SN photosphere pro- 
jection on the sky plane is an ellipse, in general a non-cancelled 
Unear polarization is expected perpendicular or paraUel to the 
major axis (Kasen et al. 2003). Therefore, if the expansion el- 
lipticity evolves with time, keeping the direction of both axes 
unchanged, then the Stokes parametres should also move on a 
straight line of the Stokes plane. Once the Milky Way Galactic 



interstellar polarization (GIP) has been corrected, this straight 
line is expected to be displaced from the origin of the Stokes 
plane due to the host galaxy interstellar polarization (HGIP). 
Using several different methods Maund et al. (2009) estimated 
a HGIP ofP~ 1.2% and ~ 134° for SN2008D/XT080109. 
The contribution of additional polarization components (variable 
line polarization, irruption of the stellar core, variable HGIP, 
etc.) can distort this rough ellipsoidal stellar picture, resulting 
in a more complex geometrical configuration (see Maund et al. 
2009, and references therein). Asphericities in explosive stellar 
deaths tend to be low, usually yielding optical polarization val- 
ues below ~ 1% (Wang & Wheeler 2008). However, counterex- 
amples exist, such as the Type Ic SN 1997X (Wang et al. 2001) 
or XRF060218/SN2006aj (Gorosabel et al. 2006; Maund et al. 
2007). 

In this paper we report optical polarimetric observations of 
SN2008D/XT080109 and SN2007uy aimed at obtaining geo- 
metrical information on the expansion of these events. In addi- 
tion we show the results of millimetre (mm) observations which 
were used to infer information about the host galaxy extinction 
and consequently about the HGIP impact on the observed polar- 
ization. 



2. Observations and data analysis 

2.1. Optical observations 

Table 1 displays the log of our optical observations. All the data 
were taken in the V-band. The observations were carried out 
with the 2.2m Calar Alto telescope (CAHA), the 2.5m Nordic 
Optical Telescope (NOT) and the 8.2m Unit Telescope 2 of the 
Very Large Telescope (VLT). 

The NOT observations were performed with ALFOSC, 
through a calcite and a 1/2 wave retarder plate. Four images of 
the field were acquired, rotating the retarder plate at 0°, 22.5°, 
45.0° and 67.5°. The calcite plates reduce the ALFOSC field 
of view (FoV) to 140" in diameter (pixel scale of O'.'19/pix). 
SN 2007uy was left out of the NOT images due to this reduced 
FoV. 

The CAHA observations were based on the CAFOS instru- 
ment. The CAFOS polarization unit uses a WoUaston prism in- 
stead of a calcite plate and has a strip mask on the focal plane 
to avoid accidental overlapping on the CCD. The total FoV of 
CAFOS is composed by 14 strips of 9' x 18" each, providing 
a large enough FoV to image both SN2008D/XT080109 and 
SN 2007uy simultaneously. 

The VLT observations were done with FORSl, with a setup 
similar to CAFOS. As in CAFOS, flie large FoV of FORSl al- 
lowed us to simultaneously image SN2008D/XT 080109 and 
SN 2007uy. In some of the FORS 1 visits we used the 1x1 bin 
read-out mode in order to avoid saturation of SN 2007uy while 
keeping a high signal to noise ratio for SN2008D/XT 080109. 
We note that, in addition to the CCD biiming, the seeing and 
transparency conditions varied significantly from night to night. 
This made it impossible to have a bright unsaturated field star 
common to all our images. Apart from SN 2008D/XT 080109, 
SN 2007uy was the object most frequently imaged, being present 
in file FoV of aU the FORSl and CAFOS images. This fact 
provided us with the opportunity of using SN 2007uy as a po- 
larimetric reference star on the Stokes plane (see discussion of 
Sect. 3.2). 
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Date" Telescope Exposure Seeing 
2008, UT (+Instrument) time (s) (//) 


SN2008D/XT 080109 
P±ap (%) S ± (Tg (deg) V ± (fy 


SN2007uy'' 
P±(Tp (%) B±as (deg) Y ± a% 


Jan 13.1644 NOT(+ALFOSC) 12x600 0.7 


0.95 ±0.20 114.9 ±5.9 18.44 ±0.05 




Tan 1 ^ 964.9 VT T^+FOR S 1 'y'\^(Y' 1 

Jail IJ.^KJ^iL V 1 V_ T^X^ WiVkJ i ^ Z,JJ\J i.U 


n(S8 + 99 11Q9 + Q'S 18'^Q + OO'S 

(0.65 ± 0.33)'' (124.5 ± 13.3) 


1 4.S + n 9^ 178 h + 4 (S IS 78 + OS 

(1.59 + 0.04) (177.4 + 0.7) 




+ 98 1 0fi 1 + 4. 1 8 98 + 0*1 




Jan 17.2518 VLT(+F0RS1) 8x285 1.0 


0.84 ±0.45 138.3 ±8.1 17.90 ±0.05 
n 14 + 51^ C141 + 72^ 


1.43 ±0.14 182.5 ±5.1 15.69 ±0.05 

n ^9 + 04^1 n774 + 071 


Tan^0 18S8 VT Tr+FOR S 1 4x98^ OS 


1 OS + OfS 1 ^ + 1 7 1 7 ^ 1 + OS 

(1.11 ±0.08) (134.0 ±2.0) 


1 + OfS 178^ + 11 1 + OS 

l.U'T- ^ U.UU 1 /O.J X i.i 1U.J>' ^ U.UJ 

(1.59 ±0.04) (177.4 ±0.7) 


Tan 30 2022 YTT^+FORSI^ 4x28S OS 


128 + 006 132 5+14 1727 + 005 
(1.21 ±0.08) (133.3 ±1.9) 


1 56 + 06 1 76 1 + 1 1 1 6 36 + 05 
(1.59 ± 0.04) (177.4 ±0.7) 


Feb 11.9326 CAHA(+CAFOS) 13x400 2.0 


1.58 ±0.97 126.6 ±16.4 17.82 ± 0.09 
(1.19 ±1.51) (131.6 + 20.5) 


0.92 ±0.88 179.8 ±26.0 17.09 ±0.11 
(1.59 ± 0.04) (177.4 ±0.7) 


Feb 26.1107 CAHA(+CAFOS) 16x600 1.9 


< 2.76' 18.52 + 0.09 
(< 2.940 


1.81 ± 0.63 173.5 ±14.4 17.56 + 0.08 
(1.59 ± 0.04) (177.4 ±0.7) 


Feb 26.9444 NOT(+ALFOSC) 16x400 1.3 


0.98 ±0.80 124.6 ±23.2 18.79 ±0.08 




Mar 2. 1102 VLT(+FORSl) 4x525 1.1 


1.42 ±0.46 139.0 ±9.1 18.80 ±0.05 
(1.17 ±0.51) (134.7 ±11.4) 


1.83 ±0.19 174.6 ±3.2 17.55 ± 0.05 
(1.59 ± 0.04) (177.4 ±0.7) 


Mar 28.0207 VLT(+FORSl) 4x525 1.1 


0.75 ±0.57 111.2 ±21.9 19.42 ± 0.05 
(0.73 ±0.78) (112.7 + 21.0) 


1.58 ± 0.42 179.6 ± 7.6 17.98 ± 0.05 
(1.59 ±0.04) (177.4 + 0.7) 



a Mean observing epoch. 

b For the NOT epochs SN 2007uy was out of the FoV, so the degree of linear polarization (P), the position angle (8), and the magnitude (V) could not be determined. 
c Calibration based on Malesani et al. (2009). The photometric eiTors include the zero point unceitainty. 

d The polarimetric cycle was interrupted and restarted several times. The SN 2007uy observing epoch and the exposure time are slightly different from the SN 20()8D/XT080109 ones, 



corresponding to Jan 13.2708 UT and 990 seconds, respectively. 
e The SN 2008D/XT 080109 P and e values given for VLT and CAHA in between parentheses assume a constant SN 2007uy polarization given by (P) = 1.59 ± 0.04% and {6} = 177.4 ± 0.7°. 
/ la- upper limit. 



The images were reduced using standard procedures running 
under IRAF' . The master flat field image was created combin- 
ing sky flat field frames taken without polarization units in the 
Ught pafli. This process was applied separately for the 3 instru- 
ments. In principle, the host galaxy background might pose a 
potential problem for the photometric accuracy of SN 2007uy 
and SN2008D/XT080109. Inspection of the NGC2770 region 
around each SN shows that SN 2007uy is surrounded by a more 
inhomogeneous background than SN2008D/XT 080109. Thus 
we double-checked the impact of the host galaxy background 
level determination by performing aperture photometry with 
radii ranging from 0.5 to 2.5 times the FWHM (Full Width at 
Half Maximum), and varying the iimer radii of the background 
annuli from 3 to 5 times the FWHM. The annuli widths were 
also varied from 2.5 to 4 times the FWHM. We verified that the 
resulting Stokes parametres were consistent within error bars, 
independent of the apertures and annuli used. Then the apertures 
and annuh yielding minimum errors were adopted (typically the 
apertures and the inner radii used were the FWHM and 3-4 times 
the FWHM, respectively). In any case, in order to account for 
potentially dismissed photometric error sources, the statistical 
analysis of Sect. 3.2.1 was duplicated considering also photo- 
metric errors augmented by 20%. 

For the ALFOSC and CAFOS data the determination of the 
Stokes parametres was done by fitting the S (0) function with the 
corresponding Milky Way GIF normahzation factor (di Serego 
Alighieri 1997). The GIF normalization factor was calculated 
using Galactic field stars. 



In the case of FORS 1 the Stokes parametres were calculated 
via Fourier expansion (see expressions 4.5 and 4.6 of FORSl 
manual'^). Then the Stokes parametres of all the objects in the 
images were corrected for the FORSl spurious polarization cor- 
rection, which considers the instrumental polarization depen- 
dence across the FoV ^. Next, the V-band angle offset of the 
1/2 wave retarder plate was included (given in the Table 4.7 
of FORSl manual). AU flie objects (including both SNe) of all 
FORSl visits used the same angle offset, minimizing a potential 
source of internal scattering in the Stokes plane. Finally, the GIF 
correction was introduced by shifting the origin of the Stokes 
plane to the bary centre of the field stars''. 

Due to the reduced ALFOSC FoV only two unsaturated 
bright field stars were available for the GIP correction of the 
NOT data. Fortunately for the three NOT epochs the stars used 
for the NOT GIP correction remained the same. This assures in- 
ternal consistency of the position of the NOT Stokes parametres 
on the Stokes plane. Thus, a potential relative shift (owing to dif- 
ferent GIP corrections) of the NOT Stokes parametres with re- 
spect to each other was minimized. This fact will be of relevance 
in the further discussion since, as will be argued in Sect. 3.2.2, 
the polarimetric data might suggest the existence of a symmetry 
axis which is consistently determined by different instrumental 
setups. 

The wider FoV of FORSl and CAFOS provided a larger 
number of stars for the GIP correction than with ALFOSC. Due 
to different observing conditions (seeing, transparency and CCD 



' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 



2 Doc. N. VLT-MAN-ESO-13100-1543. Issue 82.1, date 27/02/2008. 
^ http://www.eso.org/paranal/instruments/forsl/inst/pola.html 
The FORS 1 reduction steps explained above follow the recommen- 
dations given by Dr. T. Szeifert, priv. comm. 
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binning) at VLT and C AH A the stars varied from one night to an- 
other. This was not critical for the FORSl and CAFOS images, 
since the GIP correction for both SN2008D/XT080109 and 
SN 2007uy were always calculated using the same set of stars 
for each epoch. And most importantly, given that SN 2007uy was 
very well detected in all the VLT and CAHA images, we could 
use this object (justified by a statistical analysis, see Sect. 3.2.1) 
as a secondary calibrator and keep it fixed on the Stokes plane 
(see Sect. 3.2.2). 

It is important to stress that the GIP corrections appUed to the 
different data- sets were very similar since the Galactic extinction 
towards NGC2770 is only Ay = 0.1 mag (Schlegel et al. 1998; 
Cardelli et al. 1989). So in all cases the GIF corrections applied 
were low (~ 0.2 - 0.3%), in agreement with the Galactic polar- 
ization predictions (Pmax = 2.9% X Ay = 0.29%; (Serkowski et 
al. 1975; Cardelli et al. 1989). In fact, the GIP corrections were 
always well below the typical level of our polarization detections 
(~ 1 - 1.5%). 

Verification of the photometry, calibration and the reduction 
was done by observation of polarimetric standards. The standard 
stars observed were NGC 2024-1 and Velal-95 at VLT (Fossati 
et al. 2007; Whittet et al. 1992), HD19820 at CAHA (Wolff et 
al. 1996), and HD9485 1, BD25727, HD25 1204 and G19 1B2B at 
NOT (Tumshek et al. 1990; Schmidt et al. 1992). Table 1 shows 
the inferred and P values, once they were corrected for the 

GIP and the statistical bias (multiplying P by -^1 - (crp/P)^) 
due to the fact that P is a positive quantity (Wardle & Kronberg 
1974). The SN2008D/XT080109 P and 9 values displayed in 
Table 1 for VLT and CAHA in between parentheses assume a 
constant SN2007uy polarization given by (P) = 1.59 ± 0.04% 
and (0) = 111 A ± 0.7°. 

For some observing epochs we explored maximising the 
time resolution of our polarimetric monitoring at the expense 
of enlarging the polarimetric errors. This was done, when pos- 
sible, by splitting the polarimetric data of one night in cycles 
of four images, i.e., the minimum block of images necessary 
to get a polarimetric point. Unfortunately most of the polar- 
ization detections showed modest significance levels (between 
1.1 cr and 4.7 cr), so we could not spUt the data. Only the VLT 
data-set of January 30 could be divided in two cycles of four 
images, while still keeping high-significance polarization detec- 
tions (above 17.5 cr for each of the two cycles). We verified 
that the joint data acquired on January 30 (composed of two cy- 
cles of four images) yields Stokes parametres consistent with 
the ones obtained when the data of that night are split in two 
halves. This separation provided an extra high-quality data-point 
for the statistical analyses carried out in Sect. 3.2.1, Sect. 3.2.2 
and Sect. 3.2.3. 

The V-band polarimetry values synthesized by Maund et al. 
(2009) for SN 2008D/XT 0801 09 on January 3 1 .22 and February 
15.18 UT bring us the opportunity to cross-check our results. 
Our two closest observations, carried out on January 30.2022 
and February 11.9326 UT, yield P and values fully consis- 
tent with Maund et al. (2009). We are aware that this compara- 
tive exercise, although satisfactory, is only limited to two epochs 
close to the light curve maximum and exclusively focused on 
SN2008D/XT080109, since Maund et al. (2009) did not ob- 
serve SN 2007uy. 

2.2. Millimetre observations 

SN2008D/XT080109 was observed at 1.2 mm on January 25, 
28, and 30, 2008 with the MAMBO 11 bolometer on the IRAM 



Table 2. Log of the millimetre observations. 





Date" 


Telescope Field 


Integration 


X 


Flux density 




2008, UT 




time (min) 


mm 


[mjy/beam] 




Jan 25.931 


IRAM 30m SN2008D 


3x20 


1.2 


2.46 ± 0.54 




Jan 28.200 


IRAM 30m SN2008D 


1x20 


1.2 


2.29 ± 1.20 




Jan 30.171 


IRAM 30m SN2008D 


2x20 


1.2 


1.44 ±0.60 


Jan 25.931-30.171 IRAM 30m SN 2008D 


6x20 


1.2 


2.03 ± 0.38* 




Jan 31.925 


IRAM 30m Adjl'^ 


2x 16 


1.2 


2.39 ± 0.70 




Jan 31.936 


IRAM 30m Adj2'' 


2x 16 


1.2 


2.81 ± 0.71 




Jan 23.854 


PdB SN2008D 


1 x60 


3.3 


0.65 ±0.15 




Nov 9.336 


PdB SN2008D 


Ix 111 


2.9 


0.03 ± 0.10 


a 


Mean observin 


g epoch. 








h 


Weighted average of the three previous table Unes. 






c 


R.A.=09'' 09'" 


3a085, Dec.=+33° 08' 29'.'98 (J2000). 






d 


R.A.=09'' 09"" 


3«211, Dec.=+33° 08' 10'.'25 (J2000). 







30m Telescope (see Table 2). The data were reduced follow- 
ing the standard procedure with the mopsic data reduction soft- 
ware. In order to estimate the NGC 2770 contribution at 1 .2 mm 
to our previous SN2008D/XT080109 observations we mea- 
sured, on January 31, the flux of two adjacent regions (Adjl 
and Adj2). Such regions bracket the SN2008D/XT 080109 po- 
sition in the radial direction towards the NGC 2770 nucleus, 
botii at a distance of ~ 10" from SN2008D/XT080109 (beam 
FWHM at 1.2 mm ~ 11"). They were selected sufficiently 
close to SN2008D/XT080109 to have a good estimate of its 
background emission, but suflficiently far to avoid most of 
the SN2008D/XT 080109 flux wifliin flie beam. The adjacent 

1 .2 mm pointings yielded fluxes consistent with those obtained 
at the SN2008D/XT080109 position (see Table 2), implying 
that the 1 .2 mm data are very Ukely dominated by the NGC 2770 
background flux and not by SN 2008D/XT 080109. 

In addition, a 4.3cr detection was achieved on January 23 at 

3.3 mm with the Plateau de Bure (PdB) interferometer centered 
on SN 2008D/XT080109 (see Table 2) with an observing beam 
size of 2'.'75 X 1'.' 11. The last miUimetre observation was carried 
out on November 9.336 UT with PdB at 2.9 mm, yielding a 3cr 
upper limit of 0.30 mJy/beam. This non-detection suggests that 
the flux detected at 3.3 mm on January 23 mostly originated from 
SN 2008D/XT 080109 and not from the host galaxy dust. 



3. Results 

3. 1. Properties of the polarizing ISM towards 
SN2008D/XT080109 

The Galactic reddening towards NGC 2770 is very low (Schlegel 
et al. 1998, £(5 -V) = 0.032), implying a V-band extinction of 
only Ay = 0.1 mag (Cardelli et al. 1989). Therefore the con- 
tribution of the Galactic dust to the measured 1.2 mm flux is 
negligible. Thus, assuming that most of the 1.2 mm emission 
comes from optically thin dust emission in the host galaxy inter- 
stellar medium (ISM) towards SN2008D/XT080109, it is pos- 
sible to roughly estimate the host optical extinction as: Ay - 
(1.0865i.2)/(Bi.2(rdust) n„h)y-(Kv/Ki.2).Bi,2(Tiust) is the Planck 
function of the dust at a temperature Tdust, ^mh is the main- 
beam solid angle, 5 1.2 is the flux density per beam at 1.2 mm, 
and kv/ki.2 is the ratio of the visual extinction coefficient to 
the 1.2 mm dust opacity, which was estimated to be (kv/ki 2) = 
(4 + 2) X 10'* (Kramer et al. 1998). 
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In order to estimate Ay a value of T^ust has to be assumed. 
Domingue et al. (1999) inferred Tdust - 21 ± 2K for the 
colder dust component of a set of three spiral galaxies simi- 
lar to NGC2770. On the other hand, the spectral energy dis- 
tribution of NGC 2770 integrated over the entire galaxy yields 
Tdus, = 30 + 5 K (Thone et al. 2009). So we adopted both 
Tdust = 21 + 2 K and Tdu,t = 30 + 5 K. 

Using 5 1.2 = (2.03 ± 0.38)mJy/beam (line four of Table 2; 
weighted average of the first three table lines) and both Tdust - 
21 + 2K and Tdust = 30 + 5K, we obtain Ay = 0.43 + 0.31 
and Ay = 0.27 + 0.15, respectively. However, as 5 1.2 corre- 
sponds to the flux integrated in the beam along the projected 
thickness of the host galaxy (which includes contribution from 
any ISM located beyond the XT), we can only provide an upper 
limit to Ay. To be conservative, we will assume an extinction of 
Ay < 0.43 -H 0.31 = 0.74, obtained with Tdust = 21 + 2 K. This Ay 
value is ~ 0.5 - 1 .8 mag lower than the line-of-sight Ay values in- 
ferred by other authors based on a broad diversity of techniques 
(MazzaH et al. 2008; Soderberg et al. 2008; Malesani et al. 2009; 
Modjaz et al. 2009). An unrealistic temperature of Tdust ~ 10 K 
would be necessary to have an Ay value in agreement with the 
above authors. 

This apparent disagreement can be explained by the impos- 
sibility to resolve a clumpy host galaxy ISM with our 1.2 mm 
beam (FWHM=H"). The clumpy ISM would show three main 
properties; /) an ISM composed of cells superimposed on a lower 
extinction background, //) cells displaying typical angular sizes 
smaller than the beam size, ///) typical angular separation be- 
tween cells small enough to allow the beam to contain several 
cells. Given that the extinction determined by the 1.2 mm data 
is an average of the integrated flux received in the beam, this 
scenario would explain the higher extinction derived along the 
XT line-of-sight (i.e., spectroscopy). On the other hand, condi- 
tion ///) would explain the fact that the pointings adjacent to the 
SN2008D/XT080109 position yielded similar 1.2 mm fluxes, 
and hence extinction values. 

We can estimate a rough upper limit of the typical cell sizes 
as follows. First, we assume that the flux differences (or internal 
dispersion) between our three 1.2 mm pointings (see Table 2; 
weighted flux average of SN 2008D/XT080109, Adj 1 and Adj2) 
were entirely due to the statistical fluctuations in the average 
number of cells (A^) contained in each beam. In our case the 
fliree 1.2 mm pointings yield 2.03 + 0.38, 2.39 + 0.70 and 
2.81 + 0.71 mJy/beam, so the dispersion is 0.41 mJy/beam. If 
the cells were identical and randomly distributed, then the ratio 
between the average flux and the dispersion of the 3 pointings 
would be approximately V^V. In our case the average flux of 
the 3 pointings (calculated by weighting with the corresponding 
flux errors) is 2.24 mJy/beam, so this yields ~ 30 cells/beam. 
However, the dispersion is very likely not only due to statisti- 
cal fluctuations of A^, i.e, the dispersion probably also has an 
instrumental/calibration/photon-noise component. So the statis- 
tical fluctuations due to would be likely lower than 0.41 
mJy/beam and therefore we can only set a lower limit ofN> 30. 

Additionally, in order to match the Ay derived from the 1.2 
mm flux with the spectroscopic Ay values reported in the lit- 
erature, the emitting region (the total area covered by the cells 
within the beam solid angle) should fill only^ 10-30% of Q„,i. 
This would explain the low Ay < 0.74 mag derived using the 
1.2 mm flux (Kramer et al. 2008) for cells having a line-of- 
sight spectroscopic extinction of Ay = 1.2 - 2.5 mags. Thus, as- 
suming identical circular cells, a NGC 2770 distance of 27 Mpc 

' For dust temperatures ranging from T = 21 to T = 30K. 
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Fig. 1. Upper panel: Stokes parametres of SN 2008D/XT 080109 and 
SN2007uy, once the GIP correction has been included. The cloud of 
measurements representing SN 2007uy seems to be more compact than 
the SN2008D/XT 080109 one. The different lines represent the lin- 
ear fits to the SN2008D/XT 080109 Stokes parametres obtained when 
the CAHA (dot-dashed), VLT (continuous), NOT (dotted) and all the 
data (VLT-l-CAHA-l-NOT, long-dashed) are considered. The directions 
of all the straight lines are statistically consistent (see Table 4), which 
might indicate the presence of a preferential axis. This possible evi- 
dence is strengthened when SN2007uy is used as reference (see be- 
low). Bottom panel: Stokes parametres of SN2008D/XT 080109 and 
SN2007uy, if SN2007uy is kept fixed at its barycentre. The direc- 
tions of all the straight lines are statistically more consistent than in 
the upper panel (compare the internal consistency of the different ori- 
entations in Table 3 and Table 4), reinforcing the possible presence 
of an instrumental-independent dominant symmetry axis present in the 
SN 2008D/XT 080109 data. 



(Soderberg et al. 2008), and a Q.„,h filling factor of 30%, we can 
impose a maximum cell diameter of D < 0.16 kpc. Adopting 
a beam filling factor of 10%, we would obtain smaller cells, so 
we consider 0.16 kpc as a robust upper limit. As a comparison it 
is interesting to mention that most (~ 98%) of the giant molec- 
ular clouds resolved by millimetric observations in the nearby 
spiral M33 show sizes between 0.05 and 0.16 kpc (Rosolowsky 
et al. 2007), compatible with our observations. Similarly, most 
of the HII regions in nearby spiral galaxies show physical scales 
consistent with < 0.16 kpc (Gonzalez-Delgado & Perez 1997). 
We note that the dominant P is the result of all the material in- 
tegrated along the line-of-sight, so we can not conclude that the 
dominant P is exclusively due to the XT circumstellar dust. 
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Fig. 2. Polarimetric and photometric evolution of SN 2008D/XT 080109 (panels a,c,e) and SN 2007uy (panels b,d,f). For consistency we used the 
same symbols as in Fig. 1. The long-dashed lines indicate the position of the barycentre ((Q), ((/)) of each SN. The two bottom panels (e and f) 
show the (/-band light curves of SN2008D/XT 080109 and SN2007uy uncorrected for Galactic reddening. None of the two SNe seem to obey 
any obvious smooth long-term polarimetric evolution correlated with their optical light curves. 



3.2. Properties of SN2008D/XT080109 and SN2007uy on 
the Stokes plane 

The P and 6 values displayed in Table 1 were corrected for the 
GIP and the polarimetric statistical bias (Wardle & Kronberg 
1974). The further polarimetric analysis will deal with the prop- 
erties of Q and U of both objects on the Stokes plane, so the P 
bias correction factor (included in Table 1) will not be consid- 
ered. 

3.2.1. Is SN2008D/XT080109 polarimetrically variable? 

The Stokes parametres of SN2008D/XT080109 and 
SN2007uy, corrected for the GIP, are plotted in Fig. 1. Both 
events show polarization levels of 1 - 1.5%, significantly lower 
than the one measured in XRF060218/SN2006aj (Gorosabel 
et al. 2006; Maund et al. 2007). Assuming both a GIP law and 
a Galactic extinction law (Serkowski et al. 1975; Cardelli et al. 
1989) for NGC2770, the Ay values reported in the hterature 
(~ 1 .2 - 2.5 mag) would yield a conservative HGIP polarization 
upper limit of P,„,„ = 2.9% x Ay ~ 3.5 - 7.3%, consistent 
with the polarization measured for SN2008D/XT080109 and 
SN2007uy. 

As seen in the upper panel of Fig. 1, the cloud of points 
corresponding to SN2007uy seems to be more clustered than 
the one of SN 2008D/XT 080109. For the two stellar explosions 
we quantified the probability that the Stokes parametre distri- 



butions are consistent with no time evolution. First, we double- 
checked using Monte Carlo methods that, for the average num- 
ber of counts of our images, both Q and U follow Gaussian dis- 
tributions. This fact assures that the;^^^ test is an appropriate sta- 
tistical tool. Second, we determined the barycentre of the two 
objects on the Stokes plane {{Q}, (U); see the long-dashed lines 
of Fig. 2). Then we calculated for the Q and U of the two events 
the;^'^/dof value with respect to {Q} and {U}, respectively^. 

For SN2007uy the distribution of Q {U) shows 
;i^Vdof=4.18/7 C^f2/dof=4.32/7) with respect to {Q} ({U}). 
For SN2008D/XT080109 the distribution of Q (U) yields 
;^-Vdof=19.06/10 C^Vdof =22.96/ 10), clearly larger than the 
;^'^/dof value of SN2007uy. These values of;i'^/dof were used 
to obtain the corresponding probabilities. For SN2007uy the 
probability that Q and U are constant are 0.758 and 0.742, 
respectively. SN 2008D/XT080109 shows lower probabilities, 
0.040 and 0.011, respectively. Thus the probabilities that Q and 
U are simultaneously constant are 0.563 for SN2007uy and 
4.3 X lO-'* for SN2008D/XT 080109. 

We are aware that different eiTor sources affecting the mea- 
surements (background determination uncertainties, flat field 
coiTection inaccuracies, no instrumental polarization correction 
available along the FoV for the CAHA and NOT images, etc.) 
might still not be included in the error bars, so the values of 

' Xl/dof = ZU^ f I (k- 1) and^/dof = I^U^ f I ik- 
l),k being the number of visits and dof=fc - 1 the degrees of freedom. 
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X^/dof could be overestimated for both sources (and hence the 
above probabilities underestimated). Thus, we repeated the ex- 
ercise augmenting the photometric error bars of both objects and 
all the field stars by 20%. If we do that, the probability that both 
Q and U are constant are 0.791 for SN2007uy and 0.021 for 
SN2008D/XT080109, respectively. 

We have checked, if due to some uncontrolled reason, the 
(Q,U) distribution on the Stokes plane is just the result of mixing 
data acquired with different telescopes/instruments. In order to 
inspect the existence of this potential instrumental artifact, we 
redid the statistical analysis separating the images that contain 
both objects according to the telescope used (hence the NOT is 
excluded from this comparative study between SN 2007uy and 
SN2008D/XT080109 since SN2007uy was not imaged by the 
NOT). Thus we determined separately for the VLT and CAHA 
the probabilities of having constant Stokes parametres. 

Using only the VLT data points, we determined that the prob- 
ability that both Q and U are constant are 0.295 for SN2007uy 
and 0.007 for SN 2008D/XT080109, respectively. If we en- 
large the photometric errors by 20% the probabilities would be 
then 0.532 and 0.056, for SN 2007uy and SN 2008D/XT080109, 
respectively. Considering exclusively the CAHA data-set we 
can not reach strong conclusions given the reduced number of 
visits (2) and large error bars. However, even with the lim- 
ited CAHA data, SN 2007uy shows a higher probability than 
SN2008D/XT080109 of being polarimetrically constant, 0.457 
versus 0.403. If the CAHA photometric errors are augmented by 
20%, then the probabilities are 0.530 and 0.423 for SN2007uy 
and SN2008D/XT080109, respectively. 

Independent of the data subset considered and the additional 
photometric uncertainty introduced, SN2008D/XT080109 re- 
iteratively displays lower probabilities than SN 2007uy of be- 
ing polarimetrically constant. In fact, it is not obvious to ex- 
plain how non-accounted photometric/calibration/instrumental 
uncertainties could systematically affect only the values derived 
for SN2008D/XT080109 and not those of SN2007uy, which 
consistently shows a higher degree of clustering on the Stokes 
plane. It is interesting to note that this is against the a priori 
expectations, given that the more inhomogeneous host back- 
ground of SN2007uy (in comparison to SN2008D/XT 080109) 
might represent a potential uncontrolled source of noise in its 
Stokes parametres. All the above arguments seem to suggest that 
SN 2008D/XT080109 could show an intrinsic variable polariza- 
tion component added to the dominant HGIP. This suggestion is 
reinforced when SN 2007uy is used as a secondary polarimetric 
standard (see Sect. 3.2.2). However, it is important to note that a 
constant polarization signal from SN 2007uy is atypical of previ- 
ously observed core-collapse SNe (Wang & Wheeler 2008). We 
point out that none of the two SNe seem to obey any obvious 
smooth long-term polarimetric evolution correlated to their light 
curves (see Fig. 2). 

3.2.2. A symmetry axis for SN2008D/XT080109? 

Given that both SN2007uy and SN2008D/XT080109 were 
imaged simultaneously with VLT and CAHA, we can use 
SN2007uy as a reference star with constant Stokes parametres. 
In other words, we can assume, justified by the above probabil- 
ities, a constant projected geometry for SN 2007uy plus a con- 
stant HGIP. As an additional argument, it is also interesting to 
note that the agreement between the synthetic broad-band points 
reported by Maund et al. (2009) and our data is even better when 
SN 2007uy is fixed on the Stokes plane. 
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Fig. 3. Upper panel: The evolution of the linear fit^-/dof as a func- 
tion of the fitted line orientation when constant polarization is assumed 
for SN 2007uy. The minima correspond to the optimal orientation of the 
straight lines displayed in the Fig. 1 bottom panel, also given in Table 3. 
As seen the global shape of the evolution for the three independent data- 
sets (VLT, NOT and CAHA) are quite equivalent, their minima being 
systematically located at similar negative angles. The long-dashed line 
shows the ^-/dof evolution when all the data are jointly considered. 
We note that in the case of CAHA, given that dof=0, ;^^^/dof actually 
represents x' ■ Bottom panel: The evolution of the normalized probabil- 
ity associated with the above panel. As displayed the three independent 
data-sets show their probability maxima at a comparable orientation an- 
gle. The maximum is reinforced when all the data are combined (long- 
dashed line). Given that dof=0 for CAHA, the corresponding probabil- 
ity distribution was obtained through a Monte Carlo method. We refer 
the reader to the main text for further information. 



Thus, using SN2007uy as calibrator, a linear fit to the 
SN2008D/XT080109 data yields a satisfactory ;i^Vdof= 1.9/9 
(long-dashed straight line of Fig. 1 bottom panel). Hence, the 
SN2008D/XT080109 Stokes parametres might suggest the ex- 
istence of a dominant symmetry axis, offset from the origin due 
to the HGIP. Furthermore, as seen in the first three lines of 
Table 3 the inferred symmetry axes show a preference to neg- 
ative orientations independent of the telescope-i-instrument em- 
ployed in the observations''. This effect can be seen in the bot- 
tom panel of Fig. 1, which shows the linear fits carried out both 
when the data are considered separately (CAHA, VLT, NOT) 
and jointly (CAHA -H VLT H-NOT). The linear fits were obtained 
minimizing x^, which is defined as the weighted distance per- 
pendicular to the fitted line^ We considered this x^ definition 
because the variables to be fitted {Q,U) and their correspond- 
ing eiTors are treated symmetrically (Boogs et al. 1990; Babu & 
Feigelson 1996). 

In order to quantify the degree of consistency that the 
symmetry axes of the 3 data-sets could show (when constant 
SN 2007uy Stokes parametres are assumed) we mapped the evo- 
lution of the linear fit;i'^/dof when the orientation of the sym- 
metry axis is varied from -90 to -1-90 degrees. The upper panel 

' We define the orientation as the angle of the fitted straight line with 
respect to the horizontal axis representing Q on the Stokes plane. The 
orientation ranges between -90 and -f 90 degrees. 

** Also known as weighted total least squares (WTLS), weighted rig- 
orous least squares, or weighted orthogonal regression method (see 
Lemmerling & van HufFel 2002, and references therein). 
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Table 3. Linear fits on the Stokes plane of SN2008D/XT080109. The linear fits are shown in the bottom panel of Fig. 1, where we used SN2007uy 
as a reference star with constant Stokes parametres. Column 1 displays the telescope used and column 2 the number of visits or points on the Stokes 
plane. Column 3 provides the;i'^/dof of the optimal linear fit. Column 4 shows the orientation angles at the minima (maxima) of the Fig. 3 upper 
(bottom) panel. Columns 5 and 6 provide the lower and upper error bars around the orientation angle of the linear fit. As seen the three independent 
data-sets (CAHA, NOT and VLT) show orientations well consistent within errors. Last table line reports the results when all the data are combined. 



Telescope 


Number of 


Minimum 


Orientation angle 


Lower Error Bar 


Upper Error Bar 




visits 




(degrees) 


(degrees) 


(degrees) 


CAHA 


2 


0" 


-76.1 


49.4* 


48.5* 


NOT 


3 


0.001/1 


-64.8 


46.1 


46.7 


VLT 


6 


1.198/4 


-54.8 


52.2 


51.8 


CAHA+VLT+NOT 


11 


1.922/9 


-42.6 


23.3 


22.9 



a In this particular case, = ^ind dof={) for the minimum, so^^/dof could also be formally represented by 0/0. 
b Errors obtained through Monte Carlo methods. 



of Fig. 3 shows this evolution. We warn the reader that in the 
particular case of CAHA dof=0, so what is plotted for CAHA 
in the Fig. 3 upper panel actually represents x^- As seen all the 
three independent data-sets show their ;f^/dof minima approx- 
imately at the same orientation. The position of all the ;f^/dof 
minima were positively confirmed using the f itexy routine im- 
plemented by Press et al. (1992). The lower panel ofFig. 3 shows 
the normalized probabilities corresponding to the;k'^/dof values 
of the above panel. 

For the case of CAHA, given that dof=0, the plotted proba- 
bility evolution was obtained by simulating mock CAHA Stokes 
parametres created by a weighted Monte Carlo method assuming 
Gaussian errors for Q and U. Then given two simulated (Q, U) 
CAHA pairs the orientation of the line passing through those 
Stokes parametres is determined. The repetition of this process 
allowed us to construct the histogram of the line orientations for 
the CAHA data. Thus, the lower panel of Fig. 3 displays for 
CAHA the histogram of orientations for 10^ Monte Carlo simu- 
lations^. 

For each data- set the upper and lower error bars in the ori- 
entation angle (displayed in Table 3, colunms 5 and 6) were cal- 
culated by integrating the ±34.145% percentile of the probabil- 
ity distribution around the corresponding maximum. As seen in 
Table 3 the axes inferred from the three independent data-sets are 
consistent within errors. Hence, the inferred symmetry axes are 
statistically consistent in all cases, regardless of the instrumental 
setup used to determine it. It is also interesting to remark that 
the NOT Stokes parametres, which did not use SN2007uy as 
a secondary calibrator (so they might apply a slightly different 
GIP correction), are also satisfactorily fitted by a straight line. 
Furthermore, the orientation of the NOT linear fit agrees with 
the ones inferred based on the VLT and CAHA data. All these 
arguments strengthen the reality of this possible symmetry axis. 

The existence of a synmietry axis on the Stokes plane could 
be explained by an axisymmetric explosion where the direction 
of the symmetry axis is constant, but the eccentricity evolves 
with time. We stress that our polarimetric data can provide in- 
formation about the projected geometry on the plane of the sky 
exclusively. Thus, we can not discard that the constant Stokes 
parametres of SN 2007uy could be caused by a variable eccen- 
tricity expansion when it is viewed pole-on. 



' This Monte Carlo method, used for determining the probability dis- 
tribution for CAHA, is not optimal for the VLT and NOT data. For the 

VLT and the NOT a set of simulated Stokes parametres does not deter- 
mine a line passing through them (they contain more than 2 visits) and 
a formal fit is necessary. So, the histogram of orientations would mix 
linear fits with different qualities. 



3.2.3. The robustness of the possible SN2008D/XT080109 
symmetry axis: must we assume a constant geometry 
forSN2007uy? 

One might wonder if the possible symmetry axis of 
SN2008D/XT 080109 is just an artifact due to the assumption 
of a constant geometry for SN 2007uy. Thus we repeated the ex- 
ercise explained in Sect. 3.2.2 using the original unrecalibrated 
SN2008D/XT080109 Stokes parametres. The results are dis- 
played in the upper panel of Fig. 1 and in Table 4. Note that 
the NOT Stokes parametres, could not be recalibrated, so the fit 
is identical to the one displayed in Table 3. 

The SN2008D/XT080109 entire data set 
(CAHA+ VLT-hNOT) is stiU satisfactorily fitted (x^/dof^ 1 .066) 
by a straight line. The derived orientation i-49.S'^\lj-l deg), 
is fully consistent with the orientation determined based on 
the recalibrated Stokes parametres (-42.6;^23 3 deg). On the 
other hand, the three separated data sets yield orientations 
with internal consistency within Icr, so as in Sect. 3.2.2 the 
inferred orientations seem to be instrumentally independent. It 
is also noticeable that the recalibrated Stokes parametres yield 
orientations (-76.1, -64.8, and -54.8 deg) with lower internal 
dispersion than the unrecalibrated ones (88.9,-64.8,-86.7 deg). 

The fit quality to the VLT Stokes parametres is clearly re- 
duced with respect to the recalibrated case, but it is still fairly 
satisfactory (x^/dof - 1.375). It is worthwhile to note that the 
major source of ihisx^/dof increment comes from the lack of 
overlap within la of the SN2008D/XT 080109 Stokes parame- 
tres corresponding to VLT3a and VLT3b (red and green circles 
of the Fig. 1 upper panel), which appear slightly separated (by 
|Aei=0.0012 and |Af/|=0.0022, corresponding to a relative shift 
of 1.4cr and 2.5cr, respectively). 

One could be tempted to speculate on a possible rapid (Q,U) 
variabihty in visits VLT3a and VLT3b of SN 2008D/XT 080109. 
However, firstly the polarimetric variation between VLT3a and 
VLT3b is not statistically significant, and secondly (and most 
importantly), this internal shift between VLT3a and VLT3b is 
essentially mimicked in the same visits by the behaviour of 
SN2007uy. This suggests that this systematic shift present in 
both SNe does not reflect a real polarimetric variation, but is 
probably just a spurious polarization fluctuation due to sUghtly 
different GIP corrections for VLT3a and VLT3b which affects 
both SNe in a systematic way. In fact, in both visits the seven 
stars (identical for both visits) used for the GIP are fainter than 
SN 2007uy and SN 2008D/XT 080109, so they show larger pho- 
tometric statistical fluctuations. Thus, when using SN 2007uy as 
a calibrator the Stokes shift between VLT3a and VLT3b almost 
disappears, overlapping within Icr for SN 2008D/XT080109 
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(and obviously for SN2007uy). This example illustrates the 
importance of using the brightest common object in the field 
(SN 2007uy, even brighter than SN 2008D/XT080109) as a ref- 
erence star in order to remove internal spurious shifts below 
~ 0.2% from the Stokes plane. 

We conclude that, at a lower significance level, the unre- 
calibrated data still allow the existence of a symmetry axis for 
SN2008D/XT 080109 to be plausible. Therefore, in order to 
suggest a symmetry axis it is not strictly necessary to assume 
that the SN 2007uy geometry is constant. However, this assump- 
tion strengthens its possible existence. 

3.3. Final remarks 

Maund et al. (2009) report a variation in the position angle of the 
He I /15876 line based on the comparison of their two observ- 
ing epochs. Thus, we have to be cautious in explaining the ori- 
gin of the possible SN 2008D/XT080109 polarization variabil- 
ity. The SN2008D/XT080109 polarization variation suggested 
by the analysis of the broad-band data could not be explained, to- 
tally or partially, just by an axisymmetric expansion. We can not 
discard that the polarization variability could be, to some degree, 
caused by the changing relative strengths in the line/continuum 
emission, only one of which may be polarized, as discussed by 
Maund et al. (2009). Thus V-band polarimetric data alone can 
not resolve these potential components, and in principle only 
spectropolarimetry could study the impact of this, and other pos- 
sible effects, on the evolution of the Stokes parametres. However, 
the limited time-coverage of the spectropolarimeric data (only 
two epochs available, covering 2 weeks around the light curve 
maximum) does not make it possible to derive strong conclu- 
sions about potential eff'ects which could impact our eleven visits 
spanning ~ 75 days. We can only conclude that our data suggest 
a preference to an approximately axisymmetric configuration, 
although we can not exclude deviations from axisymmetry to- 
wards a more complex geometry, as discussed by Maund et al. 
(2009). 

Since SN2007uy occurred earlier than 
SN2008D/XT 080109, we have checked to what degree a 
possible early fast-evolving epoch of SN 2007uy may have been 
missed by our observations. An inspection of the light curves 
of the two events shows that our data cover the rising phase, 
the maximum and the decay of both objects. The polarization 
inferred for SN2007uy in these 3 phases is consistent with 
being constant. Obviously, we can not discard that SN 2007uy 
showed polarization variability earlier than our first epoch, later 
than the last epoch or in the unavoidable temporal gaps between 
observations. So, we stress that our conclusions on SN 2007uy 
and SN2008D/XT 080109 only refer to the time covered by our 
observations. 

Under the mentioned limitations, the distribution of 
SN 2007uy on the Stokes plane imposes an upper limit on the 
variability of its polarization. First, we calculated the + 1 sigma 
confidence variability interval of Q (U) around (Q) {{U}) us- 
ing the weighted variance"^ of the Q (U) sample as an statistical 
estimator. This estimator is suitable because both Q and U fol- 

As defined by the N.I.S.T. (National Institute of 
Standards and Technology; Gaithersburg, Maryland, USA, 

http://www.nist.gov/index.html), 5^ = I.loWQ,(Qi- {Q}f, 

being wq. = ai ^ . ctq- the standard deviation of g,, and N the 

number of SN2007uy visits. The expression of Su is analogous, 
replacing Qi with t/,, (Q) with (U), and ctq^ with cry;. 



low Gaussian distributions (as opposed to P). Then, the (Q,U) 
weighted variances were propagated yielding a ±34.145% con- 
fidence interval of A/'±34.i45% =!oo9^'' around the SN2007uy 
barycentre. Finally, using the analytic expression (Al) given in 
the Appendix of Wardle & Kronberg (1974) we estimated a 
±45% confidence interval of AP±4s% =^q}5%. Thus, we esti- 
mate that over the course of our observations of SN 2007uy its 
total polarization signal does not change by more than 0.29% 
(90% confidence interval). We double-checked that Monte Carlo 
methods impose similar limits on the SN2007uy polarization 
variability {AP±45% ~ 0.3%). 



4. Conclusions 

For both stellar explosions our optical polarization data seem 
to be dominated by the HGIP. This conclusion is supported 
by our 1.2 mm observations, performed at and around the 
SN2008D/XT080109 position on NGC2770. The 1.2mm 
measurements are consistent with no intrinsic emission from 
SN2008D/XT080109 at this wavelength, and can be explained 
by the host galaxy dust emission. The Ay inferred from our 

1.2 mm observations is ~ 1.2 mag lower than the line-of-sight 
extinction values reported in the literature. This apparent contra- 
diction can be solved if the NGC2770 projected dust distribu- 
tion around SN2008D/XT 080109 is composed of dense clumps 
(Ay ~ 1.2 - 2.5 mag) with a typical size < 0.16 kpc and a low 
filling factor (10-30%) in our 1.2 mm beam. 

As a bonus, we also report the detection of 0.65 mJy/beam at 

3.3 mm coincident with SN 2008D/XT080109 14 days after the 
X-ray outburst. This emission agrees fairly well with the light 
curve modeled at 95 GHz (Soderberg et al. 2008). Furthermore 
an observation at 2.9 mm carried out on November 9.3 UT 2008 
(304.8 days after the XT), imposed a 3cr flux upper limit of 0.30 
mJy/beam, so we conclude that the 3.3 mm detection is mostly 
intrinsic emission from SN2008D/XT080109. 

SN2008D/XT 080109 seemed to be surrounded by an exter- 
nal thick He mantle (MazzaU et al. 2008; Modjaz et al. 2009; 
Soderberg et al. 2008; Tanaka et al. 2009b), in contrast to flie 
He-poor XRF060218/SN2006aj (Pian et al. 2006; Soderberg et 
al. 2006; Sollerman et al. 2006). Given that the asphericity of 
core-collapse supernovae is more pronounced in the inner layers 
(Wang & Wheeler 2008), this could explain to some degree the 
modest SN 2008D/XT 080109 intrinsic polarization in compari- 
son toXRF060218/SN2006aj (Gorosabel et al. 2006; Maund et 
al. 2007). 

Notwithstanding the important HGIP, a statistical analysis 
of the distribution of flie SN 2008D/XT 080109 Stokes parame- 
tres suggests that SN2008D/XT 080109 could exhibit an intrin- 
sic variable polarization component embedded in the dominant 
HGIP. In contrast, the SN 2007uy polarization data can be de- 
scribed by a constant projected geometry. This conclusion is 
achieved independent of the data being considered jointly or sep- 
arately according to the telescope used. 

Assuming that the SN 2007uy projected geometry is constant 
in time, then the evolution of flie SN2008D/XT080109 Stokes 
parametres could show a symmetry axis. The same conclusion is 
achieved if no constant polarization is assumed for SN 2007uy, 
but at a lower significance level. This possible symmetry axis is 
also independent of the instrumentation used to perform the ob- 
servations, supporting its reality. This potential axis could be ap- 
proximately explained by an axisymmetric aspherical expansion 
with variable eccentricity. However, our broad-band polarimet- 
ric data can not exclude more complex geometric configurations. 
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Table 4. Linear fits on the Stokes plane of SN2008D/XT 080109 when no constant Stokes parametres are assumed for SN2007uy. The corre- 
sponding straight lines are shown in the upper panel of Fig. 1. Column 1 displays the telescope used and column 2 the number of visits or points 
on the Stokes plane. Column 3 provides the;^f-/dof of the optimal linear fit. Column 4 shows the orientation angles at the minima (maxima) of 
the Fig. 3 upper (bottom) panel. Columns 5 and 6 provide the lower and upper error bars around the orientation angle of the linear fit. As seen the 
three independent data-sets (CAHA, NOT and VLT) show more scattered orientations (and also larger ;k'^/dof values) than the ones displayed in 
Table 3, but they are still consistent within errors. Last table line reports the results when all the data are combined. 



Telescope 


Number of 


Minimum 


Orientation angle 


Lower Error Bar 


Upper Error Bar 




visits 


X^dof 


(degrees) 


(degrees) 


(degrees) 


CAHA 


2 


(T 


88.9 


53.9* 


48.8* 


NOT 


3 


0.001/1 


-64.8 


46.1 


46.7 


VLT 


6 


5.501/4 


-86.7 


27.9 


26.8 


CAHA+VLT+NOT 


11 


9.599/9 


-49.8 


17.4 


17.2 



a In this particular case,;tr^ = and dof=0 for tlie mimimim, so^'^/dof could also be formally represented by 0/0. 
b Errors obtained through Monte Carlo methods. 



like possible deviations from axisymmetry. Our results are con- 
sistent with the SN 2008D/XT080109 asphericity inferred from 
spectroscopic data (Modjaz et al. 2009; Tanaka et al. 2009a). We 
suggest that at least the projected, if not the intrinsic, geometry 
of the two explosive events could be different. 
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